EXACT SOLUTIONS TO THE 
NONLINEAR SCHRODINGER EQUATION 



Tuncay Aktosun and Theresa Busse 
Department of Mathematics 
University of Texas at Arhngton 
Arhngton, TX 76019-0408, USA 

Francesco Demontis and Cornehs van der Mee 
Dipartimento di Matematica e Informatica 
Universita di Caghari 

Viale Merello 92 
09123 Caghari, Italy 

Abstract: A review of a recent method is presented to construct certain exact solutions 
to the focusing nonlinear Schrodinger equation on the line with a cubic nonlinearity. With 
motivation by the inverse scattering transform and help from the state-space method, an 
explicit formula is obtained to express such exact solutions in a compact form in terms 
of a matrix triplet and by using matrix exponentials. Such solutions consist of multisoli- 
tons with any multiplicities, are analytic on the entire xt-plane, decay exponentially as 
X — i> ±oo at each fixed t, and can alternatively be written explicitly as algebraic combina- 
tions of exponential, trigonometric, and polynomial functions of the spatial and temporal 
coordinates x and t. Various equivalent forms of the matrix triplet are presented yielding 
the same exact solution. 

Mathematics Subject Classification (2000): 37K15 35Q51 35Q55 

Keywords: Nonlinear Schrodinger equation, exact solutions, explicit solutions, focusing 

NLS equation, NLS equation with cubic nonlinearity, inverse scattering transform 



1 



1. INTRODUCTION 

Our goal in this paper is to review and further elaborate on a recent method [3,4] to 
construct certain exact solutions to the focusing nonlinear Schrodinger (NLS) equation 

iut + Uxx + '^\u\'^u = 0, (1) 

with a cubic nonlinearity, where the subscripts denote the corresponding partial derivatives. 

The NLS equation has important applications in various areas such as wave propaga- 
tion in nonlinear media [15], surface waves on deep waters [14], and signal propagation in 
optical fibers [9-11]. It was the second nonlinear partial differential equation (PDE) whose 
initial value problem was discovered [15] to be solvable via the inverse scattering trans- 
form (1ST) method. Recall that the 1ST method associates (1) with the Zakharov-Shabat 
system 

dip{X, X, t) 
dx 

where u{x,t) appears as a potential and an asterisk is used for complex conjugation. 
By exploiting the one-to-one correspondence between the potential u{x, t) and the corre- 
sponding scattering data for (2), that method amounts to determining the time evolution 
u{x,0) >—>■ u{x,t) in (1) with the help of solutions to the direct and inverse scattering 
problems for (2). 

Even though we are motivated by the 1ST method, our goal is not to solve the initial 
value problem for (1). Our aim is rather to construct certain exact solutions to (1) with 
the help of a matrix triplet and by using matrix exponentials. Such exact solutions turn 
out to be multisolitons with any multiplicities. Dealing with even a single soliton with 
multiplicities has not been an easy task in other methods; for example, the exact solution 
example presented in [15] for a one-soliton solution with a double pole, which is obtained 
by coalescing two distinct poles into one, contains a typographical error, as pointed out in 
[13]. 
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—iX u{x,t) 
-u{x,t)* iX 



ip{X,x,t), 



(2) 



In constructing our solutions we make use of the state-space method [6] from the 
control theory. Our solutions are uniquely constructed via the explicit formula (9), which 
uses as input three (complex) constant matrices A, B, C, where A has size pxp, B has size 
p X 1, and C has size 1 xp, with p as any positive integer. We will refer to {A, B, C) a triplet 
of size p. There is no loss of generality in using a triplet yielding a minimal representation 
[3,4,6], and we will only consider such triplets. As seen from the explicit formula (9), our 
solutions are well defined as long as the matrix F{x^ t) defined in (8) is invertible. It turns 
out that F{xjt) is invertible if and only if two conditions are met on the eigenvalues of the 
constant matrix A; namely, none of the eigenvalues of A are purely imaginary and that no 
two eigenvalues of A are symmetrically located with respect to the imaginary axis. Our 
solutions given by (9) are globally analytic on the entire xt-plane and decay exponentially 
as a; ^ ±00 for each fixed t e R as long as those two conditions on the eigenvalues of A 
are satisfied. 

In our method [3,4] we are motivated by using the 1ST with rational scattering data. 
For this purpose we exploit the state-space method [6]; namely, we use a matrix triplet 
{A, B, C) of an appropriate size in order to represent a rational function vanishing at 
infinity in the complex plane. Recall that any rational function R{X) in the complex plane 
that vanishes at infinity has a matrix realization in terms of a matrix triplet {A,B,C) as 

R{X) = -iC{XI - iA)B, (3) 

where / denotes the identity matrix. The smallest integer p in the size of the triplet yields 
a minimal realization for R{X) in (3). A minimal realization is unique up to a similarity 
transformation. The poles of R{X) coincide with the eigenvalues of (iA). 

The use of a matrix realization in the 1ST method allows us to establish the separa- 
bility of the kernel of a related Marchenko integral equation [1,2,4,12] by expressing that 
kernel in terms of a matrix exponential. We then solve that Marchenko integral equation 
algebraically and observe that our procedure leads to exact solutions to the NLS equa- 
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tion even when the input to the Marchenko equation does not necessarily come from any- 
scattering data. We refer the reader to [3,4] for details. 

The explicit formula (9) provides a compact and concise way to express our exact so- 
lutions. If such solutions are desired to be expressed in terms of exponential, trigonometric 
(sine and cosine), and polynomial functions of x and t, this can also be done explicitly 
and easily by "unpacking" matrix exponentials in (9). If the size p in the matrices A, B, 
C is larger than 3, such expressions become long; however, we can still explicitly evaluate 
them for any matrix size p either by hand or by using a symbolic software package such 
as Mathematica. The power of our method is that we can produce exact solutions via (9) 
for any positive integer p. In some other available methods, exact solutions are usually 
tried to be produced directly in terms of elementary functions without using matrix expo- 
nentials, and hence any concrete examples that can be produced by such other methods 
will be relatively simple and we cannot expect those other methods to produce our exact 
solutions when p is large. 

Our method is generalizable to obtain similar explicit formulas for exact solutions 
to other integrable nonlinear PDEs where the 1ST involves the use of a Marchenko inte- 
gral equation [1,2,4,12]. For example, a similar method has been used [5] for the half-line 
Korteweg-de Vries equation, and it can be applied to other equations such as the modified 
Korteweg-de Vries equation and the sine-Gordon equation. Our method is also general- 
izable to the matrix versions of such integrable nonlinear PDEs. For instance, a similar 
method has been applied in the third author's Ph.D. thesis [8] to the matrix NLS equation 
in the focusing case with a cubic nonlinearity. 

Our method also easily handles nonsimple bound-state poles and the time evolution of 
the corresponding bound-state norming constants. In the literature, nonsimple bound-state 
poles are usually avoided due to mathematical complications. We refer the reader to [13], 
where nonsimple bound-state poles were investigated and complications were encountered. 
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A systematic treatment of nonsimple bound states has recently been given in the second 
author's Ph.D. thesis [7]. 

Our main results are summarized in Section 2 and some explicit examples are provided 
in Section 3. For the proofs, further results, details, and a summary of other methods to 
solve the NLS equation exactly, we refer the reader to [3,4]. 

2. MAIN RESULTS 

In this section we summarize our method to construct certain exact solutions to the 
NLS equation in terms of a given triplet {A, B, C) of size p. For the details of our method we 
refer the reader to [3,4]. Without any loss of generality, we assume that our starting triplet 
{A^B^C) corresponds to a minimal realization in (3). Let us use a dagger to denote the 
matrix adjoint (complex conjugate and matrix transpose), and let the set {a^j^j^ consist 
of the distinct eigenvalues of A, where the algebraic multiplicity of each eigenvalue may be 
greater than one and we use rij to denote that multiplicity. We only impose the restrictions 
that no aj is purely imaginary and that no two distinct Oj values are located symmetrically 
with respect to the imaginary axis on the complex plane. Let us set A-,- := iaj so that we 
can equivalently state our restrictions as that no will be real and no two distinct Aj 
values will be complex conjugates of each other. Our method uses the following steps: 

(i) First construct the constant p x p matrices Q and N that are the unique solutions, 
respectively, to the Lyapunov equations 



QA + A^Q = C^C, 



(4) 



AN + N A^ = BB\ 



(5) 



In fact, Q and can be written explicitly in terms of the triplet [A, B, C) as 
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N^ — [ dX{XI -iA)-^BB\XI + iA^)-\ (7) 

where 7 is any positively oriented simple closed contour enclosing all in such a way 
that all A* lie outside 7. The existence and uniqueness of the solutions to (4) and (5) 
are assured by the fact that Xj ^ X* for all j = 1, 2, . . . , m and Xj 7^ A^ for k ^ j. 

(ii) Construct the p x p matrix valued function F{x,t) as 

F{x, t) := e2^'--«(^')'* + Q e-^^^-^^^'^N. (8) 

(iii) Construct the scalar function u{x, t) via 

u{x, t) := -2B^F{x, t)-^CK (9) 

Note that u{x, t) is uniquely constructed from the triplet {A, B, C). As seen from (9), 
the quantity u{x,t) exists at any point on the xt-plane as long as the matrix F{x,t) 
is invertible. It turns out that F{x,t) is invertible on the entire xt-plane as long as 
Xj ^ X* for all j = 1, 2, . . . , m and A^- 7^ A^ for k ^ j. 

Let us note that the matrices Q and N given in (6) and (7) are known in control 
theory as the observability Gramian and the controllability Gramian, respectively, and 
that it is well known in control theory that (6) and (7) satisfy (4) and (5), respectively. In 
the context of system theory, the invertibility of Q and N is described as the observability 
and the controllability, respectively. In our case, both Q and N are invertible due to the 
appropriate restrictions imposed on the triplet {A, B, C), which we will see in Theorem 1 
below. 

Our main results are summarized in the following theorems. For the proofs we refer 
the reader to [3,4]. Although the results presented in Theorem 1 follow from the results 
in the subsequent theorems, we state Theorem 1 independently to clearly illustrate the 
validity of our exact solutions to the NLS equation. 
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Theorem 1. Consider any triplet {A,B,C) of size p, corresponding to a minimal repre- 
sentation in (3), and assume that none of the eigenvalues of A are purely imaginary and 
that no two eigenvalues of A are symmetrically located with respect to the imaginary axis. 
Then: 

(i) The Lyapunov equations (4) and (5) are uniquely solvable, and their solutions are 
given by (6) and (7), respectively. 

(ii) The constant matrices Q and N given in (6) and (7), respectively, are selfadjoint; i.e. 

= Q and = N. Furthermore, both Q and N are invertible. 

(iii) The matrix F{x, t) defined in (8) is invertible on the entire xt-plane, and the function 
u{x,t) defined in (9) is a solution to the NLS equation everywhere on the xt-plane. 
Moreover, u{x,t) is analytic on the entire xt-plane and it decays exponentially as 
X —>■ ±oo at each fixed t e R. 

We will say that two triplets (A, B, C) and (A, B, C) arc equivalent if they yield the 
same potential u{x,t) through (9). The following result shows that, as far as constructing 
solutions via (9) is concerned, there is no loss of generality is choosing our starting triplet 
(A, S, C) of size p so that it corresponds to a minimal representation in (3) and that all 
eigenvalues aj of the matrix A have positive real parts. 

Theorem 2. Consider any triplet {A, B, C) of size p, corresponding to a minimal repre- 
sentation in (3), and assume that none of the eigenvalues of A are purely imaginary and 
that no two eigenvalues of A are symmetrically located with respect to the imaginary axis. 
Then, there exists an equivalent triplet {A, B, C) of size p, corresponding to a minimal 
representation in (3) in such a way that all eigenvalues of A have positive real parts. 

The next two result given in Theorems 3 and 4 show some of the advantages of using a 
triplet (A, S, C) where all eigenvalues of A have positive real parts. Concerning Theorem 
2, we remark that the triplet {A, B, C) can be obtained from [A, B, C) and vice versa with 
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Q 



POO 

/ ds [Ce 



the help of Theorem 5 or Theorem 6 given below. 

Theorem 3. Consider any triplet {A,B,C) of size p, corresponding to a minimal repre- 
sentation in (3). Assume that all eigenvalues of A have positive real parts. Then: 

(i) The solutions Q and N to (4) and (5), respectively, can be expressed in terms of the 
triplet {A,B,C) as 

poo 

-^^]t[Ce-^^], N= ds[e-^'B][e-'^'B]'^. (10) 

Jo 

(ii) Q and N are invertihle, selfadjoint matrices. 

(iii) Any square submatrix ofQ containing the (l,l)-entry or {p,p)-entry ofQ is invertible. 
Similarly, any square submatrix of N containing the (l,l)-entry or {p,p) -entry of N 
is invertible. 

Theorem 4. Consider a triplet {A, B, C) of size p, corresponding to a minimal repre- 
sentation in (3) and that all eigenvalues aj of the matrix A have positive real parts and 
that the multiplicity of aj is rij. Then, there exists an equivalent triplet {A,B,C) of size 
p corresponding to a minimal representation in (3) in such a way that A is in a Jordan 
canonical form with each Jordan block containing a distinct eigenvalue aj and having —1 in 
the superdiagonal entries, and the entries of B consist of zeros and ones. More specifically, 
we have 





r^i 


. 


- 











A2 . 









A = 








, B = 






. 


. 


-^m - 




-Bm- 



(11) 
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Aj:= 








aj .. 
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, Bj := 
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.. 






.1. 



C,:=[c 



Cji CjoJ, 



where Aj has size rij xuj, Bj has size rijXl, Cj has size Ixrij, and the (complex) constant 
Cj(^nj-i) is nonzero. 



We will refer to the specific form of the triplet {A,B,C) given in (11) as a standard 
form. 

The transformation between two equivalent triplets can be obtained with the help of 
the following two theorems. First, in Theorem 5 below we consider the transformation 
where all eigenvalues of A are reflected with respect to the imaginary axis. Then, in The- 
orem 6 we consider transformations where only some of the eigenvalues of A are reflected 
with respect to the imaginary axis. 

Theorem 5. Assume that the triplet {A,B,C) of size p corresponds to a minimal realiza- 
tion in (3) and that all eigenvalues of A have positive real parts. Consider the transfor- 
mation 

{A, B, C, Q, N, F) ^ {A, B, C, Q, N, F), (12) 

where {Q,N) corresponds to the unique solution to the Lyapunov system in (4) and (5), 
the quantity F is as in (8), 

A=-A\ B = -N-^B, C=-CQ-\ Q = -Q-^, N = -N-\ 

and F and u are as in (8) and (9), respectively, hut by using {A, B,C,Q, N) instead of 
{A, B,C,Q, N) on the right hand sides. We then have the following: 

(i) The matrices Q and N are selfadjoint and invertible. They satisfy the respective 
Lyapunov equations 

f QA + A^Q = C^C, 

\ ^ _ _ _ (13) 

[ AN-\-NA^ = BB^. 

(ii) The quantity F is transformed as F = Q~^FN~^ . The matrix F is invertible at every 
point on the xt-plane. 

To consider the case where only some of eigenvalues of A are reflected with respect to 
the imaginary axis, let us again start with a triplet {A, B, C) of size p and corresponding 
to a minimal realization in (3), where the eigenvalues of A all have positive real parts. 
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Without loss of any generality, let us assume that we partition the matrices A, B, C as 

C=[Ci C2], (14) 





'A, 


" 






A = 




, B = 









A2 




B2 



so that the q x q block diagonal matrix Ai contains the eigenvalues that will remain 
unchanged and A2 contains the eigenvalues that will be reflected with respect to the 
imaginary axis on the complex plane, the submatrices Bi and Ci have sizes q x 1 and 
Ixq, respectively, and hence A2, B2, C2 have sizes {p — q)x{p — q), {p — q)xl, lx{p — q), 
respectively, for some integer q not exceeding p. Let us clarify our notational choice in 
(14) and emphasize that the partitioning in (14) is not the same partitioning used in (11). 
Using the partitioning in (14), let us write the corresponding respective solutions to (4) 
and (5) as 



Q 



'Qi 


Q2 








N2 




N = 






Qz 


Q4. 






.A^3 


N4_ 



(15) 



where Qi and Ni have sizes q x q, Q4 and A^4 have sizes {p — q) x {p — q), etc. Note that 
because of the selfadjointness of Q and N stated in Theorem 1, we have 

Q\ = Qu QI^Qs, qI^Qa, nI = Ni, iV] = iV3, ivi = iV4. 

Furthermore, from Theorem 3 it follows that Qi, Q4, Ni, and A^4 are all invertible. 

Theorem 6. Assume that the triplet {A,B,C) partitioned as in (I4) corresponds to a 

minimal realization in (3) and that all eigenvalues of A have positive real parts. Consider 
the transformation (12) with {A,B,C) having similar block representations as in (14), 
{Q,N) as in (15) corresponding to the unique solution to the Lyapunov system in (4) and 

(5), 

ii = Ai, A2 = -Al Bi = B^-N2N^^B2, B2^-N^^B2, 
Ci = Ci — C2Q4 ^Qsi C2 — —C2Q4 ^, 

and {Q,N) given as 



Qi = Qi — Q2Q4 Qzi Qi = —Q2Q4 , Qs = —Q4 Qs, Q4 = —Q4 1 
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and F and u as in (8) and (9), respectively, hut by using {A, B,C,Q, N) instead of 
{A, B, C, Q, N) on the right hand sides. We then have the following: 

(i) The matrices Q and N are selfadjoint and invertible. They satisfy the respective 
Lyapunov equations in (13). 



(ii) The quantity F is transformed according to 



F = 







/ 



and the matrix F is invertible at every point on the xt-plane. 
(iii) The triplets {A,B,C) and {A,B,C) are equivalent; i.e. u{x,t) = u{x,t). 

3. EXAMPLES 

In this section we illustrate our method of constructing exact solutions to the NLS 
equation with some concrete examples. 

Example 1. The well-known "n-soliton" solution to the NLS equation is obtained by- 
choosing the triplet {A, B, C) as 



^ = diag{ai, a2, . . . ,an}, ^^ = [1 1 ... 1], C = 



Ci C2 ... Cn\, 



where aj are (complex) nonzero constants with positive real parts, B contains n entries, 
and the quantities Cj are complex constants. Note that diag is used to denote the diagonal 
matrix. In this case, using (8) and (10) we evaluate the (j, /c)-entries of the n x n matrix- 
valued functions Q, N, and F{x,t) as 



Njk = 



1 



aj + al ' 



Qjk — 



a* + ak 



-2a,,x— 4ia„t 



{a* + as){as + al)' 
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where 5jk denotes the Kronecker delta. Having obtained Q, N, and F{x,t), we construct 
the solution u{x,t) to the NLS equation via (9) or equivalently as the ratio of two deter- 
minants as 

u{x, t) 





Ct F{x,t) 



(16) 



det F{x,t) 

For example, when n = 1, from (16) we obtain the single soliton solution 

-8ct(Re[ai])2e-2«i^+4^(«i)'* 

^ ~ 4(Re[ai])2 + |2 g-4a:(Re[ai])+8t(Im[a2]) ' 

where Re and Im denote the real and imaginary parts, respectively. From (1) we see that if 
u{x, t) is a solution to (1), so is e^^u{x, t) for any real constant 9. Hence, the constant phase 
factor e*^ can always be omitted from the solution to (1) without any loss of generality. 
As a result, we can write the single soliton solution also in the form 

u{x,t) = 2Re[ai] e-^^^^^I^il+^^^^^'t^^Jsech (2Re[ai](a; - 4tlm[ai]) - log |ci/2Re[ai]|) , 

where it is seen that u{x,t) has amplitude 2Re[ai] and moves with velocity 4Im[ai]. 
Example 2. For the triplet {A, B, C) given by 





"2 " 




"l" 












-1 




1 



C^\l -1 



we evaluate Q and N explicitly by solving (4) and (5), respectively, as 



(17) 



N = 



1/4 1 
1 -1/2 

and obtain F{x,t) by using (8) as 

_ + 1 

F{x,t) = 



, Q = 



1/4 -1 
-1 -1/2 



p4x-16it _ p2x-4it i J_ -4a;-16it ^ -ix-lUt , \ 2x-Ut 

^16 4 ^2 



1 ^-Ax-16it _ 1 ^x-Ait 
4 2 



g-2a;-4it _ ^-Ax-l&it _|_ _ ^2x-m 

4 



Finally, using (9), we obtain the corresponding solution to the NLS equation as 

8e^»*(9e-^^ + 16e^^) - 32e^^»*(4e-2^ + 96^^) 



u{x^ t) = 



-128 cos(12t) + 4e-6a^ + 166^2' + 81e-2a^ + 64e2a= " 
12 



(18) 



It can independently be verified that u{x,t) given in (18) satisfies the NLS equation on 
the entire xt-plane. 

With the help of the results stated in Section 2, we can determine triplets {A, B, C) 
that are equivalent to the triplet in (17). The following triplets all yield the same u{x,t) 
given in (18): 



(i) A = 



2 
1 



B 



9/ai 
-4/^2 



C — [ai 0.2], where cki and 0:2 are arbitrary (com- 
plex) nonzero parameters. Note that both eigenvalues of A are positive, whereas only 
one of the eigenvalues of A in (17) is positive. 



(ii) A = 



-2 
1 



, B = 



16/(9a3) 
-4/(9^4) 



, C = [as q;4 ] , where and a4 are arbitrary 
(complex) nonzero parameters. Note that the eigenvalues of A in this triplet are 
negatives of the eigenvalues of A given in (17). 



(in) A = 



2 
-1 



, B = 



1/^5 
-1/q!6 



, C = [a^ chq] , where and aQ are arbitrary 



(complex) nonzero parameters. Note that A here agrees with A in (17). 



(iv) A = 









-1 



B = 



IG/cty 
-9/^8 



C = [ CKy as] , where ay and ag are arbitrary 
(complex) nonzero parameters. Note that both eigenvalues of A are negative. 

(v) Equivalent to (17) we also have the triplet {A, B, C) given by 



"9 

(1 - aQ){aQ - 2) 

QJlO 



"10 
3 — CKg 



C = [ail 



"12 



B = 



5a^Qaii + Q;ioQ!i2 ~ ^aQaiQai2 
14aioaii — Sagctioaii + 10q;i2 — lbaQai2 + 5agQ;i2 



'^lo'^ii + 3q:ioQ:iiq:i2 — 2q!9q;ioQ!iiQ!i2 + 2q!^2 ~ 3q;9Q!^2 + '^9'^i2 ' 
where ag, . . . , 0:12 are arbitrary parameters with the restriction that aiQaiiai2 7^ 0, 
which guarantees that the denominator of B is nonzero; when a^ = we must have 
0(iiOii2 7^ and choose cug as 2 or 1. In fact, the minimality of the triplet {A,B,C) 
guarantees that B is well defined. For example, the triplet is not minimal if aiiai2 = 0. 
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We note that the eigenvalues of A are 2 and 1 and that A here is similar to the matrix 
A in the equivalent triplet given in (i). 

Other triplets equivalent to (17) can be found as in (v) above, by exploiting the 
similarity for the matrix A given in (ii), (iii), and (iv), respectively, and by using (3) to 
determine the corresponding B and C in the triplet. 
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